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DUCT  LINER  OPTIMIZATION  FOR  TURBOMACHINERY 
NOISE  SOURCES 

By 

Harold  C.  Lester  and  Joe  W.  Posey 
INTRODUCTION 

In  recent  years  much  research  effort  has  been  directed  towards 
developing  mathematical  models  for  understanding  the  propagation  and 
attenuation  of  sound  in  acoustically  lined,  axi symmetric,  aircraft  engine 
ducts.  This  paper  discusses  the  related  problem  of  predicting  optimal 
multisectioned  acoustical  liners  for  circular  ducts  by  the  application  of 
numerical  optimization  alogrithms  to  analytical  duct  models. 

An  optimal  liner  is  achieved  by  judiciously  selecting  the  impedance 
characteristics  of  the  liner  so  that  maximum  sound  attenuation  is  produced 
for  a  given  design  condition.  The  often  cited  work  of  Cremer  (ref.  1) 
is,  perhaps,  the  earliest  attempt  to  define  an  optimal  liner  criterion. 
Cremer 's  approach  is  based  on  maximizing  the  attenuation  rate  of  the 
lowest-order  (least  attenuated)  acoustical  mode  propagating  in  an 
infinitely  long,  two-dimensional  duct.  Tester  (ref.  2)  recently  generalized 
and  extended  Cremer 's  results  for  application  to  rectangular  and  circular 
cross-sections  and  to  arbitrary  (higher-order)  modes.  This  extension 
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without  flow.  Rice's  model  is  based  on  the  superposition  of  a  finite 
number  of  soft-wall  radial  duct  modes.  Contour  maps  were  employed  to 
determine  optimal  liner  impedance  values. 

Apparently  Wilkinson  (ref.  4)  was  the  first  to  calculate  optimal 
liner  impedance  values  by  using  a  numerical  optimization  procedure.  He 
utilized  an  integral -equation  duct  propagation  model  coupled  with  the 
method  of  steepest  descent  (gradient  method).  His  results  for  uniform 
cylindrical  liners  with  a  plane-wave  source  compare  favorably  with  Rice 
(ref.  3).  Limited  success,  due  to  numerical  difficulties,  was  obtained 
for  a  uniform  annular  duct  liner  with  flow  and  a  bellmouth  termination. 
Quinn,  in  a  recent  paper  (ref.  5),  used  a  finite  difference  solution  of 
the  convective  wave  equation  although  details  of  his  minimization  method 
are  not  mentioned.  Results  are  presented  for  a  variety  of  axi symmetric 
multi  sectioned  liner  configurations  with  and  without  uniform  flow.  In 
another  recent  paper  (ref.  6),  Beckemeyer  and  Sawdy  investigated  the 
properties  of  optimal  multisectioned,  two-dimensional  duct  liners.  Their 
theoretical  model  is  based  on  modal  superposition  (mode  matching  method). 

A  conjugate  gradient  algorithm  was  utilized  in  the  optimization  process. 
Simple  acoustic  sources,  expressed  in  terms  of  mode  amplitude  and  phase 
angles,  were  investigated. 

In  this  paper  optimal  circular  duct  liners  are  calculated  by  using 
Zorumski's  (ref.  7)  multisectioned  theory  and  the  optimization  algorithm 
of  Davidon,  Fletcher,  and  Powell  (ref.  8).  Uniform  liners  and  liners 
with  two  and  three  sections  are  investigated.  Whereas  previously 
published  results  have  emphasized  optimization  for  plane-wave  sources,  here 
the  source  models  include  a  point  source  located  axi symmetrically  on  the 
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duct  centerline  and  a  source  model  representing  rotor/outlet-stator  viscous 
wake  interaction  effects  for  a  typical  research  compressor  configuration 
operating  at  an  axial  flow  Mach  number  of  0.4.  For  this  latter  source, 
the  classic  Kemp-Sears  theory  is  used  to  describe  the  unsteady  blade 
loadings  (ref.  9).  Optimal  liner  properties  predicted  for  the  exact  set 
of  source  amplitudes  and  phase  angles  are  compared  with  those  for 
equi parti tion-of-energy,  zero-phase  and  least-attenuated-mode  source 
approximations.  A  uniform  flow  velocity  profile  is  assumed. 

SYMBOLS 

duct  liner  radius 

axial  acoustical  energy  flux  at  z 
frequency,  m  =  Z-jif 
denotes  imaginary  part  of  {  } 
v/T 

wave  number 

circumferential  wave  number 
denotes  real  part  of  {  } 
nondimensional  admittance,  6  =  +  jBj 

angle  of  incidence 
transmission  loss  in  decibels 

DISCUSSION 

For  the  present  study,  the  mode  summation  propagation  theory  developed 
by  Zorumski  (ref.  7)  for  multi  sectioned,  axi symmetric  ducts  is  employed. 
This  approach  allows  for  the  straight  forward  modeling  of  any  of  the  liner 
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configurations  shown  in  figure  1,  which  are  made  up  of  a  sequence  of  circular 
sections  each  having  uniform  admittance  properties.  A  finite  section  of  an 
infinite  duct  or  a  truly  finite  duct  can  be  easily  handled.  Arbitrary 
source  spatial  distributions  are  defined  by  giving  the  amplitudes  and  phase 
angles  of  the  various  modes.  Also,  source  impedance  and  termination  impedance 
conditions  can  be  imposed.  All  calculations  presented  here  are  based  upon 
a  wave  source,  such  that  the  right-moving  wave  at  the  source  plane  is 
specified.  Configurations  a  and  b  (fig.  1)  are  used  for  the  point 
source  calculations,  while  configurations  c,  d,  and  e  are  used  with 
turbomachinery  sources. 

A  comparison  of  this  theory  with  an  independent  modal  theory  (ref.  3) 
applied  to  the  configuration  shown  in  part  (a)  of  figure  1,  with  a  uniform 
pressure  profile  at  z  =  0,  shows  excellent  agreement  (not  shown  here  for 
the  sake  of  brevity).  Good  agreement  also  exists  with  a  finite  difference 
solution  (ref.  5)  which  employed  a  uniform  pressure  source  at  z  =  0  and 
a  uniform  axial  impedance  at  z  =  L  equal  to  the  characteristic  impedance 
of  the  medium. 

The  optimization  process  is  indicated  in  figure  2.  For  any  given 
duct  geometry,  flow  field,  and  spatial  distribution  of  a  pure-tone  noise 
source  of  frequency  f  (with  a  time  dependence  exp  (-jZ-rrft)),  the 
Davidon-Fletcher-Powel 1  minimization  algorithm  is  used  to  predict  the  liner 
admittance  giving  the  largest  possible  transmission  loss  <j)  between  the 
axial  stations  z  =  0  and  z  =  L.  The  liner  is  taken  to  be  locally 
reacting,  so  that  its  properties  are  completely  specified  by  giving  the 
normal  acoustic  admittance,  6  =  (complex  acoustic  normal  velocity)/(complex 
acoustic  pressure),  for  each  of  its  segments.  The  transmission  loss  is 
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the  reduction  in  axial  acoustic  energy  flux  expressed  in  dB.  All  results 
which  follow  are  for  transmission  losses  across  a  length  of  the  duct  equal 
to  its  diameter  (L  =  2b). 

Figure  3  is  a  plot  of  the  maximum  transmission  loss  as  a  function  of 
normalized  frequency  kb  for  a  uniform  liner  (fig.  la)  with  a  point  source 
located  on  the  duct  axis  at  z  =  0.  Since  the  situation  is  axi symmetric 
only  circumferential  wave  number  m  =  0  modes  are  excited.  Along  with 
the  numerical  optimization  results,  a  curve  determined  from  Cremer's 
least-attenuated-mode  criterion  (ref.  2)  is  shown.  At  very  low  frequencies, 
when  only  the  lowest-order  mode  is  cut  on,  Cremer's  criterion  works  well. 
However,  significant  differences  in  the  two  predictions  occur  as  soon  as 
the  first  higher-order  radial  mode  begins  to  propagate.  Cremer's  criterion 
predicts  a  maximum  (p  at  kb  =  20  of  less  than  4  dB,  while  about  10  dB 
is  actually  possible.  This  is  simply  because  higher-order  modes  are  much 
more  strongly  excited  by  the  point  source  at  high  frequencies. 

Since  one  generally  has  difficulty  in  forming  a  clear  mental  picture 
of  acoustical  propagation  expressed  as  the  sum  of  many  modes,  it  is  of 
value  to  compare  the  high  frequency  behavior  of  the  numerically  optimized 
liner  of  figure  3  with  the  predictions  of  geometrical  acoustics.  When 
kb  »  1 ,  geometrical  acoustics  is  applicable  to  the  point  source  in  a 
short  duct  and  a  ray-tracing  solution  is  a  valid  approximation.  Therefore, 
the  acoustic  pressure  and  particle  velocity  at  any  point  on  the  duct  cross 
section  at  z  =  L  may  be  found  by  summing  the  contributions  from  the  direct 
ray  and  all  the  reflected  rays  passing  through  that  point  (see  fig.  4). 

From  this  viewpoint,  the  best  that  any  liner  could  do  would  be  to  completely 
absorb  all  incident  rays,  thereby  reducing  the  acoustic  wave  at  z  =  L 
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to  only  the  direct  radiation.  For  L  =  2b  only  12  percent  of  the  rays 
emitted  from  the  right  side  of  the  source  are  radiated  directly  through 
the  z  =  L  cross  section,  implying  an  upper  limit  for  the  transmission 
loss  of  9.6  dB  (kb  »  1).  Figure  3  shows  the  optimal  4)  predicted  from 
the  modal  theory  settling  down  to  about  this  value  as  kb  approaches  20. 

Naturally,  no  uniform,  point  reacting  liner  could  provide  complete 
absorption  of  all  incident  rays,  since  the  angle  of  incidence  0  and  the 
distance  from  the  source,  b/cos  0,  are  not  constant.  Complete  absorption 
occurs  only  when  the  surface  admittance  in  the  direction  of  the  incident 
ray  is  equal  to  the  specific  acoustic  admittance  of  a  spherical  wave  at 
the  appropriate  distance  from  the  source.  Normalized  with  respect  to  the 
characteristic  admittance  of  the  medium,  this  specific  admittance  is 

(1.0  +  j  — ).  Because  the  rays  are  not  normal  to  the  locally  reacting 

liner,  the  requried  normal  admittance  is  not  this  spherical  wave  admittance 
but  is  this  quantity  multiplied  by  cos  0  (ref.  10). 

In  estimating  the  optimal  uniform  liner  admittance  from  geometrical 
acoustics  considerations,  it  is  reasonable  to  assume  that  any  ray  at  z  =  L 
which  has  been  reflected  off  the  liner  more  than  once  makes  a  negligible 
contribution  to  the  acoustic  disturbance.  Figure  4  shows  a  ray  which  is 
reflected  only  once  and  then  just  clears  the  end  of  the  duct.  Such  a 
ray  strikes  the  liner  at  z  =  L/3  and  any  ray  which  hits  the  duct  wall 
closer  to  the  source  will  suffer  multiple  reflections.  Therefore,  using 
a  single  reflection  criterion  requires  maximizing  absorption  for  rays  which 
first  contact  the  liner  between  z  =  L/3  and  z  =  L.  For  such  rays. 
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the  average  value  of  cos  e  is  approximately  0.7,  implying  an  optimal 
normal  admittance  of  about  (0.7  +  j  ). 

The  real  and  imaginary  parts  of  the  optimal  admittance  which 
correspond  to  the  maximum  transmission  loss  of  figure  3  are  shown  in 
figures  5a  and  5b,  respectively,  as  functions  of  normalized  frequency  kb. 
These  values  determined  from  the  modal  theory  are  in  good  agreement  with 
geometrical  predictions  for  large  kb. 

Optimal  liners  with  two  sections  of  equal  length  (fig.  lb)  were 
also  calcualted  for  a  point  source.  The  optimization  procedure  was 
initialized  at  the  optimal  uniform  admittance  solutions.  The  resulting 
maximum  (j)  (fig.  6)  is  little  better  than  that  for  the  uniform  liner 
either  at  very  low  kb  values  (where  only  the  lowest-order  mode  is 
cut  on)  or  at  large  kb  (where  geometrical  acoustics  applies).  However, 
at  some  intermediate  frequencies  more  than  20  dB  improvement  is  possible. 
The  largest  segmentation  advantage  for  this  source  in  this  duct  is  in 
the  kb  range  between  4  and  6,  where  only  the  lowest  two  radial  modes 
are  cut  on  in  a  hard  walled  duct.  The  results  at  large  kb  values 
are  again  consitent  with  geometrical  acoustics,  with  the  optimal 
admittances  (not  shown)  approaching  the  single  reflection  criterion 
values  of  approximately  (0.8  +  j  ^^)  and  (0.6  +  j  ^i^).  respectively. 

When  a  flow  is  introduced  in  the  minus  z  direction,  sound 
propagation  in  the  positve  z  direction  is  retarded,  giving  a  liner 
of  fixed  length  more  time  to  attenuate  a  disturbance.  Figure  7  shows 
the  effect  of  flow  on  transmission  loss  (p  and  admittance  g  for 
upstream  propagation,  kb  =  6,  and  a  point  source.  It  can  be  seen 
that  about  twice  as  much  loss  is  possible  at  Mach  0.5  as  with  no  flow. 


7 


While  the  optimal  admittance  is  seen  to  change  considerably  with  flow, 
the  performance  of  the  zero  flow  optimal  liner  in  the  presence  of 
flow  is  near  optimal  up  to  a  Mach  number  of  about  0.3.  Thus,  for  this 
specific  case,  the  sensitivity  of  liner  performance  to  normal  admittance 
is  not  excessive. 

A  12-inch  research  compressor  currently  installed  in  an  anechoic 
chamber  at  Langley  Research  Center  is  selected  as  a  realistic  turbo¬ 
machinery  noise  source.  The  design  parameters  appropriate  to  this 
axial  compressor, configured  with  a  single  stage  having  a  rotor  with 
19  blades  and  26  outlet  guide  vanes, are  supplied  to  a  computer  program 
(ref.  9)  which  employs  the  Kemp-Sears  viscous  wake  interaction  model 
to  determine  appropriate  amplitudes  of  the  acoustic  modes  for  a 
rotor/stator  interaction  noise  source.  A  uniform  axial  flow  of 
0.4  Mach  number,  corresponding  to  a  normalized  blade  passage  frequency 
of  kb  =  16.53,  is  assumed.  According  to  the  usual  Tyler-Sofrin  analysis, 
three  m  =  -7  radial  modes  are  cut-on. 

Results  are  summarized  in  figure  8.  The  bar  labeled  "exact"  shows 
the  maximum  transmission  loss  determined  for  the  calculated  source 
structure  in  the  duct  configurations  shown  in  parts  c,  d,  and  e  of 
figure  1.  While  for  the  point  source  (as  well  as  for  the  plane-wave 
source,  though  it  is  not  discussed  here)  very  little  is  gained  from 
segmentation  at  a  kb  as  large  as  16,  this  turbomachinery  noise  can  be 
reduced  an  additional  3  dB  by  using  two  segments  and  5  dB  by  using  three 
segments . 
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The  other  bars  in  figure  8  correspond  to  variations  from  the  exact 
calcualted  source.  Optimal  attenuations  are  shown  for  one-,  two-,  and 
three- segment  liners  for  each  source  variant:  (1)  putting  equal  energy 
in  each  of  the  three  propagating  radial  modes  with  the  calculated  relative 
phasing  retained,  (2)  setting  the  amplitudes  of  the  two  higher-order 
radial  modes  to  zero,  and  (3)  replacing  the  three  complex  modal  coefficients 
with  their  absolute  values.  It  is  clear  that  while  segmentation  results 
in  a  significant  improvement  regardless  of  the  source  type,  the  liner 
cannot  be  optimally  designed  unless  both  the  modal  distribution  of 
source  energy  and  the  relative  phasing  among  these  modes  is  known. 

CONCLUDING  REMARKS 

The  results  o-^  the  present  optimization  procedure  compare  well 
with  least-attenuated  mode  theory  at  low  frequencies  and  with  geometrical 
acoustics  at  high  frequencies.  This  favorable  agreement  lends  credibility 
to  the  results  at  intermediate  frequencies,  where  there  is  no  simple  theory 
with  which  to  compare. 

It  is  found  that  the  potential  benefits  of  liner  segmentation  for 
the  attenuation  of  turbomachinery  noise  is  greater  than  would  be  predicted 
from  point  source  results.  Furthermore,  effective  liner  design  requires 
precise  knowledge  of  the  circumferential  and  radial  source  distributions. 
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ADMIHANCEa 


Figure  1.-  Duct  liner  configurations. 


LINERS 


Figure  2.-  Liner  optimization  design  concept. 


Figure  5.  -  Optimal  liner  properties  for  a  uniform  circular  liner 
(fig.  1-a)  with  a  point  source  and  L  =  2b. 


(b)  Imaginary  part  of  admittance,  R{B} 


% 


40 


TRANSMISSION 

LOSS 

4>.dB 


O 

_L. 

ai 


1.0  0- 

o 

0.8  -  o 


OPTIMAL 
ADMITTANCE  0.4 


3 


0.2 


O 


■0  2  °  ° 

I _ L 


0  0.1 


O  OPTIAAALP 
O  3=  1.05+  j0.15 

O  @  3 

_ I _ I _ \- 

0.2  0.3  0.4 

MACH  NUMBER 

O  R|3i 
□  -II3t 


□  □ 

□ 

_ I _ \ _ \ _ I 

0.2  0.3  0.4  0.5 

AAACH  NUMBER 


Figure  7.-  Optimal  properties  of  a  uniform  circular  liner  (fig. 

with  a  uniform  flow  profile,  a  point  source,  L  ■  2b 


kb  =  6. 
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and 
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Figure  8.-  Maximum  transmission  loss  of  a  research  compressor  for 
different  source  models  -  circular  liners  (fig.  1-c,  d 
kb  =  16.53,  Mach  number  =0.4  and  m  =  -7. 


